
ti.y <- "Prior"

#logRo
pr <- theta[1,1]
x <- seq(theta[1,2],theta[1,3],0.1)
if(theta[1,5]==1){
	#unif
	prdens <- dunif(x,theta[1,6],theta[1,7])
	}
if(theta[1,5]==2){
	#unif
	prdens <- dnormf(x,theta[1,6],theta[1,7])
	}

if(theta[1,5]==3){
	#Beta
	prdens <- dbeta(x,theta[1,6],theta[1,7])
	}
if(theta[1,5]==4){
	#gamma
	prdens <- dgamma(x,theta[2,6],theta[2,7])
	}

plot(x,prdens,type="h",col=3)


#Steepness
pr <- theta[2,1]
x <- seq(theta[1,2],theta[1,3],0.1)
if(theta[2,5]==1){
	#unif
	prdens <- dunif(x,theta[1,6],theta[1,7])
	}
if(theta[2,5]==3){
	#Beta
	prdens <- dbeta(x,theta[2,6],theta[2,7])
	}
